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Abstract 

Massless neutrinos are exchanged in a neutron star, leading to long 
range interactions. Many body forces of this type follow and we resum 
them. Their net contribution to the total energy is negligible as compared 
to the star mass. The stability of the star is not in danger, contrary to 
recent assertions. 
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It has been recently claimed that the multibody exchange of massless neutri- 
nos renders neutron stars unstable, as the induced self-energy exceeds the mass 
of the star [Q. Multibody neutron potentials are generated. The putative cul- 
prits would be the long distance, infrared, effects associated to neutrino exchange 
among 4 neutrons or more. Specifically, the more neutrons are involved in a given 
potential, the larger contribution to the energy would result. This would be an 
exceptional case in physics where the main contribution would not be dominated 
by the combinatorial effect of few body potentials. 

Long range forces mean infrared effects. The infrared scale of the problem 
is the radius R of the star. Additional parameters of the problem which may 
play a significant role are the strength of the weak interactions, Gp, the neutron 
density, n„, and a possible neutrino density, riy-. 

= 1.1710-^GeV-^ n„~0.4fm~^ 

~ 10 km, riy ~ 4. 10"^^ fm"^ (1) 

We first analyze the problem in the Hartree-Fock approximation: correlations 
among neutrons are neglected. That is, it will be assumed that the neutrons are 
uniformly distributed in the star. 

1 Neutrino spectrum in a neutron star 

The properties of neutrinos and antineutrinos that propagate in a medium differ 
from those in the vacuum. In particular, the vacuum energy-momentum relation 
for massless fermions, u; = |g|, where uj is the energy and |g| the magnitude of 
the momentum vector, does not hold in a medium 0] . 

Assume for the time being that the material of the neutron star is made 
exclusively of neutrons, among which neutrinos are exchanged. The density- 
dependent corrections to the neutrino self-energy result, at leading order, from 
the evaluation of Z°-exchange diagrams between the neutrino and the neutrons 
in the medium, with the propagator evaluated at zero momentum. They can 
be summarized by the following dispersion relation 

u:=\q\±h, (2) 

where 

h ~ -^/2GFnn/2 0.210"^GeV ~ -10"^ fm"^ = -10^1 A"^ (3) 

In eq. the upper (lower) sign refers to neutrinos (antineutrinos). h resums 
the zero-momentum transfer interaction of a massless neutrino with any number 
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of neutrons present in the media. Sensibly enough, it depends on the neutron 
density instead of on the total number of neutrons. 

The dispersion relation shows a displacement of the energy levels for the dif- 
ferent modes, a negative shift for neutrinos, and a positive one for antineutrinos. 
The Dirac see level is displaced. Would the neutron star occupy the whole uni- 
verse, it would just mean a change of variables, with no physical consequence. 
The finite size of the star changes the picture. Notice that b acts as the depth of 
a potential well Q It is repulsive for antineutrinos and attractive for neutrinos, 
which could condense {b plays also the role of a chemical potential). 

The effective propagator corresponding to eq.(|]) can be written as 




i= (Qo- b)-fo - g7, (5) 

an infrared safe propagator. 

The above results could be obtained as well in an effective lagrangian ap- 
proach. Since we are interested in long distance effects, and correspondingly low 
momentum exchange, it is a good approximation to use the following effective 
lagrangian, as done in ref. 



Cef[{x) = i u{x) - bwjou , (6) 

This effective lagrangian is valid inside the star. The second term summarizes 
the interactions with the neutrons, with the neutron current reduced in average 
to its charge density. For simplicity, we have assumed the neutron density to be 
constant in the star. 

The effective lagrangien is not written in terms of quarks but in terms of a 
neutron density . Its validity is thus reduced to a momentum range larger than 
a typical hadronic size, ~ Ifm. When treating long distance effects, ~ this 
implicit scale is irrelevant. 

^ When a neutrino-sea is present |^ , the above remarks hold with b given by 
b ~ tV2Gf {nn - n,,)/2 - ^^^^'^ n^ < E^, >, 

oTTi ^ 

where < -E^ > is the medium average of the neutrino energy, defined in the rest frame of the 
medium. These corrections to the fermion propagation would give higher order effects, and we 
disregard them in the present paper. Direct contributions of the neutrino sea to the energy 
level of the stars will be considered later, though. 

^ The connection with Fischbach's notation is given by the replacement b = V2Gi;-a„n„. 
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Figure 1: Diagrammatic representation of eq. (8). The simple line represents 
the neutrino propagator in real vaccum. The double line stands for the "dressed" 
neutrino propagator, eq (4), which embodies the interactions with the medium. 



Eqs. (ID and (H) are the correct solution to the neutrino field equation inside 
the star which follows from eq. Since the problematic potential energy 

under consideration corresponds to possible long distance divergences, it is worth 
to recall a well-known fact: eqs. (§) and (j^) can be derived as well in perturbation 

by summing 1,2 oo interactions of an "undressed" neutrino with the neutrons. 

In the infrared, for \q\ < b, this perturbative method fails, though, because the 
infinite sum does not converge, it is infrared divergent. The derivation of eqs. (^) 
and from eq. has no such limitation. It has incorporated the interactions 
with the media, including the long distance ones. This suggests that no long 
distance divergence is to be expected from physical effects stemming from the 
interior of the star. We could stop the argument here. For the sake of comparison 
with previous literature, we explicit the computation below. 

2 Energy density from neutrino exchange 

A traditional way to estimate the energy induced by neutrino exchange is to 
compute first the possible 2,3,4, ... exchange potentials and then add their con- 
tributions integrated over the neutron positions in the star. We use instead a 
simpler and more direct method, formally equivalent to the latter, as proven in 
section 4. 

As remarked by Fischbach, eq. (B2) in Schwinger has provided the tools 
to compute the density of weak interaction energy, w, due to neutrino exchanges. 
It is given by the difference between the energy density for a neutrino propagating 
in the "vacuum" defined by the neutron star, |0 >, and the corresponding one 
for the real, matter- free, vacuum, |0 >, 




(7) 
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where 7^o(0) is the free hamiltonian density for the propagating neutrino, and 
7i(0) includes the interaction with the neutrons, as given by eq. 

In diagramatic form, depicting the "dressed" massless fermion propagator, eq. 
(P, by a double line, the formal eq. (|^) corresponds to the computation of the 
diagrams in fig. 1, 



(2^ 



'-i)Tr 



go7o \ 



1-75 



(8) 



This is a formal expression. A regularization procedure must be chosen. Let 
us consider Pauli-Villars. The problem has no infrared divergences, and it will 
control ultraviolet ones. Upon Pauli-Villars regularization, we can make a shift 
of variables, and the total result is zero, using the symmetry —Qfj.- 






Figure 2: The diagrams in fig. developped according to Schwinger-Dyson. 

See eq. (10). 



We can gain further insight by using the Schwinger-Dyson expansion. 



1 111 

-1^ = 1 + 1 o-io -, 



(9) 



and eq.(p|) is equivalent to the computation of the diagrams in fig. 2, 



w 



J 27r 4^ ' 



(2' 



^ 1 T 



d^q 



-z)Tr 



1 T 

qolo 7^70 7 ^ 



+ 



d^q 



'-i)TT 



qolo 7^70 7^70 



(10) 



6 



The first term on the right hand side of this equation corresponds in fact 
to the self-energy of one neutron, and strictly speaking it should be substracted 
(upon Pauli-Villars regularization, it gives a null contribution by itself, though). 
We see then that the neutrino contribution to the energy density stems formally 
from the computation of the diagram b) in fig. (2), that is. 



w 



pv {2txY 



'-i)Tr 



?o7o 7^70 4&7o 7 L 



0, 



(11) 



where PV stands for Pauli-Villars. The result holds in any regularization scheme 
invariant for shifts of the energy variable. 



3 Energy density from neutrino condensate 

Before turning to the connection with Fischbach's work and the discussion of the 
regularization procedure, we should consider the presence of a neutrino sea inside 
the neutron star. This will be related to finite size effects. Indeed, Smirnov and 
Vissani have recently recalled that such a sea is present 0] , and they argued that 
the two-body potential in the star is modified ||^ due to the blocking effects of 
the sea, damping the long range forces. 

We are interested in resumming the contribution to the energy from multibody 
exchange, to all orders. In our formalism this is easily done. It suffices to replace 
in eqs. ®-(|ll]) the dressed propagator, eq.(|), by 



-^2'Kii5{f)d{-qQ)d{q,-h) 



(12) 



where the Fermi sea contains the neutrino states that have a negative energy in 
the star, i.e., according to eq. (|), |g| < |6| or 6 < go < 0. As stated in 
these neutrinos are trapped inside the star by the attractive potential while the 
antineutrinos are repelled away from the star. 

Denoting by Ws this new contribution to the energy density, it follows that 



(2^ 



1 
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'?o7o 7&7o 2TTi0{q ) 6* (-go) 6* (go - b) 670 7 L 



This equation has an inmediate transparent interpretation: it is the contribu- 
tion to the total energy density due to the neutrino Fermi-sea, which is supposed 
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to be filled up in a neutron star. Quantitatively, a tiny one. Its physical meaning 
stems from the finite size of the star, that is, from the difference of the energy 
levels between the dense medium and the free vacuum. In other words, would 
the neutron star occupy the whole universe, antineutrinos could not have been 
expelled out through its surface. Indeed, it is easy to check that Ws would vanish 
if a term representing an antineutrino sea would be added to the propagator, eq. 



Given the result in the previous section, the total energy density, eq. (|^, is 
then given by 



w = w. 



-10"^^GeVfm-l (14) 



Comparison with Fischbach's analysis 



In order to compare with Fischbach analysis eq.([Tl|) can be developed using 
the Schwinger-Dyson equation (||). This leads to 



k=2 



k 



90707(^707) ^ 



(15) 



Each term in the series of eq. (|T3p correspond exactly to those considered by 
Fischbach, with the difference that the sum does not stop at the total number of 
neutrons in the star, but runs up to 00. It should be so: a given neutrino can 
interact several times with the same neutron, a fact not included in Fischbach 
treatment. 

He has estimated the fc— body potential for non-zero momentum transfer, and 
then integrated over the relative positions of the neutrons. This is equivalent to 
computing the loop in eq. (^) in which all of the k neutron current insertions 
have vanishing momentum transfer. With the method used in this paper, an in- 
sertion of the neutrino energy operator, go7o, is present as well. Upon integration 
by parts0 with a proper regularization, it is straightforward to rewrite it as 



^ = 2^1: 77;-\i(-^)Tr 



fc=2 



kJ (27r)^ 




(16) 



The factor 1/k appearing in this formula is the one required to obtain exactly 
the expansion of the logarithm, eq. (B34), in 

Taken separately, each term in the expansion (|l^) with 4 or more neutron 
insertions. A; > 4, is infrared divergent, the degree of divergence rising with 

''This integration by parts corresponds to the steps described in eqs. (B27)-(B29) in 
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the number of neutron insertions. They would lead to a contribution to the 
energy density of the form where bR ^ I.IO^^ ^ i_ jhis reasonma; 

lead Fiscbach to conclude that neutron stars are unstable if the neutrinos are 
massless0. Nevertheless, as remarked above, the correct behaviour which follows 
from the field equations for the lagrangian, eq.(^, proves that the sum of the 
series is infrared convergent, as evident from eq. ([TI|). We are facing a situation 
where the behaviour of the different terms in a perturbative expansion differs 
essentially from the resummed, non-perturbative, result^. Multibody massless 
neutrino exchange do not destabilize the star. 

For the sake of comparison with Fischbach's paper it is worth to comment 
upon the treatment of ultraviolet divergences, even if the latter are irrelevant to 
the long distance problem under discussion, as already remarked in |T[ . As finally 
such a problem does not exist, they could be the leading contribution to the star 
self-energy, although a harmless one. 

Fischbach has correctly remarked that only terms with 4 or less neutron inser- 
tions, /c < 4, are ultraviolet divergent. He explicitly computed their finite contri- 
bution to the star self-energy. For A; = 2 he obtained a non-zero result, contrary 
to our results above. The issue is not the form of the 2-body exchange poten- 
tial, but the energy density derived from it upon spatial integration, where the 
ultraviolet divergence appears. The difference stems from the regulator choice. 
For simplicity, we have chosen above to work with Pauli-Villars, a perfectly con- 
sistent one. Fischbach implemented a cutoff non invariant for the shift of energy 
variables, a harder one, leading to a finite result. He points out that there is a 
natural ultraviolet cutoff in the problem - the "hard core" r^. 



Tc ~ 0.5 fm, (17) 

which prevents neutrons from "piling up" in space. It can be interpreted as the 
natural hadronic cutoff for the effective lagrangien, eq. ®, discussed above. To 
consider the effects of such a cutoff is tantamount to go beyond the Hartree Fock 
aproximation, where correlations among neutrons are neglected, r^ is then to be 
added to our list of physically relevant parameters in eq. ([^). Its exact shape 
is not known from first principles. We do not wish to discuss here Fischbach's 
choice, as it only affects the multibody potentials involving 4 neutrons or less, 
irrelevant for the infrared problem under discussion. The physical consequences, 
if present, should be derivable from the underlying theory, the standard model. 
We do not follow this line of research here, as the main focus of this paper are 

^Smirnov and Vissani, revised version may 24th, argue that due to Pauh blocking, the 
infrared dangerous parameter 0{bR) in Fischbach's expansion is replaced by a parameter 0(1) 
and the series might be resummed. We have performed the sum with (section 3), or without 
(section 2) Pauli blocking and find no infrared divergence. 

^From an algebraic point of view the behavior is similar to that of the function 1/(1 — x), 
for X > 1, which is wellfined even if its expansion in x diverges. 
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long distance effects, for which Vc is irrelevant. His result, even if larger than 
the contribution of the Fermi sea derived here, does not threaten the stability of 
the star. Would we use Fischbach's cutoff, our dominant contribution would be 
provided by the term with the strongest ultraviolet divergence, k = 2: from two 
body neutrino exchange. The result is, eq. (4.2) in|] 

= ^ 10-16 GeVfm-^ (18) 

87r^r^ 

to be compared to the neutron mass density ~ 0(1 GeV fm~^). Eq. (^) does 
not take into account the presence of the neutrino sea. As the latter only corrects 
the long distance behaviour, the minor changes resulting from its inclusion do not 
change the physical result for the star energy density. 



5 Conclusion 

We have shown that the resummation of multibody neutrino exchange in a neu- 
tron star results in an infrared well-behaved contribution to the star energy den- 
sity. This holds whether the presence of a neutrino sea is taken into account 
or not. The resulting energy density is negligible as compared to the star mass 
density. The star remains stable, even if neutrinos are massless. 
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